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Abstract 

We elucidate the asymptotics of the L*-quantization error induced by a sequence of L^- 
optimal n-quantizers of a probabihty distribution P on M'^ when s > r. In particular we show 
^ that under natural assumptions, the optimal rate is preserved as long as s < r + d (and for 

' every s in the case of a compactly supported distribution) . We derive some applications of these 

"j^ , results to the error bounds for quantization based quadrature formulae in numerical integration 

on M'' and on the Wiener space. 
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m 

O ! 1 Introduction 

\o 

, Optimal quantization is devoted to the best approximation in (P) (r > 0) of a random vector X : 

I (r2,^,P) M°' by random vectors taking finitely many values in R"' (endowed with a norm || . ||). 

a • When X G L''(P), this leads for every n > 1 to the following n- level L''(P)-optimal quantization 

^ . problem for the random vector X defined by 

> 



en,r{X) := inf - q{X)\l, q-.R"^ ^ ]R^,Borel function, card(g(]R^)) < n} (1.1) 



■ One shows that the above infimum can be taken over the the Borel functions 

q-.R'^ ^a:=q{R'^), acR'^, carda < n 
which are some projection following the nearest neighbour rule on their image i.e. 

{Va{a))aea being a Borel partition of R'^ satisfying 

Va{a) C IxeR'^ : ||x - a|| = min ||x - 6|| L 
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The set a = q{M.'^) is (also) called a Voronoi n-quantizer and one denotes 

X" := q{X). 

Then, if d{x, a) := minaga — a|| denote the distance of x to the set a, one has 

which shows that Cn^ri^) actually only depends on the distribution P = of X so that 



en,ri^) = en,r(^) = inf ( / d{x,aYdP{x) 

card(cj)<n \J 



The first two basic results in optimal quantization theory are the following (see 0): 

- The above infimum is in fact a maximum: there exists for every n > 1 (at least) one L^(¥)- 
optimal n-quantizer a* and if supp(P) is infinite card(a* ) = n. 

- Zador's Theorem: If L^'^''(P) i.e. J^d \\x\\'^~^'^ dP {x) < +00 for some f] > 0, then 

limn^e„,^(P) = {QriP))- e M+. 

A more explicit expression is known for the real constant Qr{P) (see 1)2. 3|) below). In particular, 
Qr{P) > if and only if P has an absolutely continuous part (with respect to the Lebesgue measure 
\d on R"*). When P has an absolutely continuous part, sequence (an)n>i of n-quantizer is U-rate 
optimal for P if 

limsupnd / d(x, a„)''P^ (dx) < +cxd 
n Jm.<i 

and is asymptotically L'" -optimal if lim^^^^^^ — ^— ^ — ^} .^^ — ^ = 1. 

Our aim in this paper is to deeply investigate the (asymptotic) L'^-quantization error induced 
by a sequence (a„)„>i of L'^-optimal n-quantizers. It follows from the monotony of s ^ || .\\s 
that {an)n>i remains an L^-rate optimal sequence as long as s < r. As soon as s > r no such 
straightforward answer is available (except for the the uniform distribution over the unit interval 
since the sequence ((^|^)i<fc<n)n>i is L'"-optimal for every r > 0). 

Our main motivation for investigating this problem comes from the recent application of optimal 
quantization in Numerical Probability to numerical integration (see ^) and the computation of 
conditional expectation {e.g. for the pricing of American options, see |1U|). Let us consider for the 
sake of simplicity the case of the error bound in the quantization based quadrature formulas for 
numerical integration. Let F : — > M be a function with a Lipschitz continuous differential 
DF. It follows from a simple Taylor expansion (see ^j) that for any random vector X with 
distribution P = ¥^ quantized by X" (a C 



E(F(X)) - E(F(X")) - E{DF{X'^).{X - AT")) < [DF]up^X - A°|2. 

If a is an L^-optimal (or quadratic) quantizer then it is stationary (see ^0] or [3]) so that 

X" = E(X|X") 
which makes the first order term vanish since 

E(L'F(X").(X - X")) = E(DF(X").E(X - X" | X")) = 0. 
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Finally, if {an)n>i is a sequence of quadratic optimal n-quantizers 

E{FiX)) - E(F(X")) < [DF]uAen,2{P)? ~ [DF]up^^. (1.2) 



nd 



Now, if the Hessian D^F does exist, is p-Holder (/9 G (0,1]) and computable, the same approach 
yields 



E(F(X)) -E(F(X"")) -E((X-X"")*i?^F(X"")(X-X"")) < - 1"" ||^+;;. (1.3) 

Consequently evaluating the asymptotic behaviour of \\X — X°'"\\2+p = (Jj^d d(x, a„)^+''(iP(x)) 2+p 
is necessary to evaluate to what extend the quadrature formula in H1.3() does improve the for- 
mer one (|1.3|) . Similar problems occur when evaluating the error in the first order quantization 
based scheme designed for the pricing of multi-asset American options or for non-linear filtering 
(see 12], J2j). One also meets such mismatch problems in infinite dimensions when dealing with 
(product) functional quantization on the Wiener space in order to price path-dependent European 
options (see the example in Section and ) . 

The paper is organized as follows: in Section [21 a lower bound for the L*(P)-quantization rate 
of convergence of an asymptotically L^'-optimal sequence (ari)n>i of n-quantizers is established. It 
shows in particular that for absolutely continuous distributions P with unbounded support, it is 
hopeless to preserve the quantization rate n~d rate in as soon as s > r + d. When s < r + d it 
provides a lower bound which can be finite. We conjecture that, when (an)n>i is L^-rate optimal 
the lower bound is in fact the sharp rate. In Section |21 we show several natural criterions on the 
distribution P which ensure that (a„)„>i is L'^-rate optimal for a given s£ {r,r + d) or even for 
any s G {r,r + d). Our criterions are applied to many parametrized families of distributions on R"^. 
We investigate by the same method in Section |3 the critical case s = r + d and the super-critical 
case s > r + d. In Section El we show that for compactly supported distributions on the real line 
the lower bound obtained in Section |21 does hold as a sharp rate. Finally, in Section |H1 we apply 
our results to the evaluation of errors in numerical integration by quantization based quadrature 
formulas in finite and infinite dimensions. 

Notations : • || • || will denote a norm on M'^ and B{x,r) will denote the closed ball centered 
at X with radius r > (with respect to this norm), d{x,A) will denote the distance between xG M.'^ 
and a subset A C M'^. 

• Arf will denote the Lebesgue measure on M*^ (equipped with its Borel cj-field B(R.'^)). 

• Let (a„)„>o and {bn)n>o be two sequences of positive real numbers. The symbol a„, x 6„ is 
for o„ = Oihn) and bn = 0{an) whereas the symbol ~ hn means an = bn + o{bn) as n ^ oo. 

• [x] is for the integral part of the real number x. 

• f (X g means that the functions / and g are proportional. 



2 The lower estimate 

In this section we derive a precise lower bound in the (r, s)-problem for non-purely singular proba- 
bility distributions P and asymptotically //''-optimal quantizers. It is expected to be best possible. 
Let r G (0,oo). Let P be a probability measure on (R"', i3(]R'^)) satisfying 

/ \\x\\''dP{x) < +00 (2.1) 
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and supp(P) is not finite. Then en,riP) £ (0,00) for every n and en,r{P) ^ as n ^ 00. A 
sequence (a„)„>i of quantizers is called asymptotically ^''-optimal for P if carda„ < n for every 
n and 

/ d{x,anYdP{x) enriPy as n ^ 00. (2.2) 

Let = f.Xd denote the absolutely continuous part of P with respect to Xd and assume that 
/ 11x11^"'''' (iP(x) < +00 for some 77 > 0. Then by the Zador Theorem (see 0) 

hm n^/^enAPy = Qr{P) (2-3) 

n — 5- 00 

where 

Qr{P) := Jr,d / /"/("+^')(iA, G [0, 00) (2.4) 

and 

Jr,d := inf n^/'^e„,,(t/([0, 1]^))^' G (0,oo). 

{U{[0, l]"^) denotes the uniform distribution on the hyper-cube [0, 1]"'). Note that the above moment 
assumption implies that J^d f'^^^'^^^^dXd < +00. 

Furthermore, for probabilities P on M'^ with P"^ 7^ 0, the empirical measures associated to an 
asymptotically L''"-optimal sequence (an)n>i of n-quantizers satisfies (see |3] Theorem 7.5 or '3 ) 

^J2 6a^Pr (2.5) 
Tl — 

where Pr denotes the L'"-point density measure of P defined by 

j!d/(d+r) 

Pr:=fr.Xd with /r := j^d/(d+r)dXd ^^'^^ 

Note that the limit Qr{P) in the Zador Theorem reads 

Qr{P) = Jr,d I f^^'^dP^ 

The quantity that naturally comes out in the (r, s)-problem, r, sG (0, 00), is 

QrAP) ■■= Js,d I fr^'^dP'^ (2.7) 

(j f/i^+^)dXd\" I /^-^/('^+^)ciA,G(0,+oo]. 
Km-i J J{f>o} 



J. 



s,d 



Theorem 1 Let r,s G (0,oo). Assume P^ ^ and [ ||x||-''dP(x) < +00 for some , > 0. Let 
(«n)n>i be an asymptotically L^ -optimal sequence of n-quantizers for P. Then 

liminfn"/'^ / d(x, a„)"dP(x) > Qrs(P)- (2.8) 

n — >oo J ' 
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Prior to the proof, let us provide a few comments on this lower bound. 

Comments. • The main corollary that can be directly derived from Theorem ^ is that 

r j,i-s/id+r)^^^ = +00 ^ lim n'/^ I d{x,anydP{x) = +00 
i{/>0} "^"^ J 

since then Qr^s{P) = +00. 

By contraposition, a necessary condition for an asymptotically L''-optimal sequence of quan- 
tizers (a„) to achieve the optimal rate n~*/'^ for the L*-quantization error is that Qr^s{P) < +00. 
But, under the moment assumption of Theorem 1 one has the following equivalence holds true 

QrAP) < +00 ^ / f-^rdP"^ = [ f^-^^dXd < +00 (2.9) 

J J{f>0} 

since J^^ f^^'^+^UXa < +00 (see ^F). 

In turn, for probability measures P satisfying \d{f > 0) = +00 a necessary condition for the 
right hand side of 1)2. 9() to be satisfied is that 

s<d + r. (2.10) 

Indeed, if s > d + r, the following chain of inequalities holds true 

d+L / \ ^ 

A,(/>0) = I l{f:^o}f-'dP'^ < [I (l{/>o}r')'^dP'') ^ = f-'/^'^'^dxA ' 

where we used that p 1-^ || . ||LP(pa) is non-decreasing since P°'{M.'^) < 1. 

On the other hand, still when s < d+r, the following criterion holds for the finiteness of Qr,s{P)'- 

3i? > 0, / ||xf^/("'+"-")+''dP(x) < +oo\ =^ Qrs{P) < +00. (2.11) 

Set p = 1 — -ipp^ G (0, 1) and u = + '&■ Then 1)2. 11(1 follows from the regular Holder inequality 

applied with p=- = 3#r^ and g = 



p d+r—s y 1— p d+r— s' 

/ fdXd < (f {f{xr\\x\rfdXd{xy ( f ||x||-"^«'dArf(x)^ ^ 

Jb{0,1)'= \JB(0,1)= / 1^5(0,1)= 



fix)\\xrdXd{x)] ([ \\x\\-^p/^^-pUXd{x)] <+oo 



using the moment assumption in (|2.1H) and up/{l — p) = d + jr" > d. 

. It is generally not true in the general setting of Theorem □ that lim n^/^ / d(., a„)^dP(x) = 

Qr,s{P) (see Counter- Example 2 in Section However, one may reasonably conjecture that this 
limiting result holds true for sequences (a„) of exactly ^''-optimal n-quantizers. Our result in one 
dimension for compactly supported distributions (see section ISJ supports this conjecture. 

• In any case, note that (|2.8() improves the obvious lower bound 

liminfn"/'^ / dix^anYdPix) > liminf .(^')' > QsiP)- 



5 



(The right inequahty is always true and comes from the proof of the Zador Theorem, see jSj). In 
fact, one even has that, for every r, sG (0, +00), 

QrAP) > Qs{P)- 

Furthermore, this inequahty is strict when r ^ s (except if / is A^-a.e. constant on {/ > 0}). Set 
p = [d + s)/s > 1, q = {d + s)/d > and a = ds/{d + r){d + s), b= {d + r- s)d/{d + r){d + s). Then 
the Holder inequality yields 



J J{f>0} 

< (^j pPdX^ (j f'^dX^ ("<" if and /^^ are not proportional) 



{/>o} 

/(d+s) ( . \d/{d+s) 



(//--..A.)' ■ 

(Q.,,(P))^. 



Proof of Theorem 1. First note that, the r-moment assumption on P implies the finiteness of 

d 

J fd+rdXd- The existence of at least one asymptotically L'"-optimal sequence («„)«>! follows from 

d ~ 

the moment assumption on P (as well as the finiteness of J^d dA^j. For every integer m > 1, 
set 

fm--= E "^''^ Er = [-<f<^]nB{0,m). 

The sequence {fm)m>i is non-decreasing and converges to /l{o</<+oo} = / Xd-a.e.. Let := {k£ 
{0, . . . , m2™ — 1} : Xd{EJP') > 0}. For every fee I^, there exists a closed set C satisfying 

Let Em £ (0, 1] be a positive real number such that the closed sets := {x E M*^ : d(x, A^) < 
}, ke Im, satisfy 

f -d_ r _d_ 

/ dAd < (1 + 1/m) / fd+^dXd < +00. 



Set 



It is clear that 



m2™--l , 
V —1 
k=0 ^ 



{fm^fm}C U (i?rW)- 

0<fc<m2'"-l 



Hence 



so that 



Xdi{fm^U})< E 



m2™-l ^ 



k=0 



32™ m? 



m>l 
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i.e., for Arf-a.e. x, fm{x) = fm{x) for large enough m so that fm converges to / Xd-a.e.. Finahy, as 
a result fm^fm^f and fm converges to / Xd-a.e.. Then, for every n > 1, 

/ {d{x,an)ydP{x) > {d{x,an)Y Jm{x)d\d{x) 



n 



^ E ^ Jd{x,an)yd\d{x). (2.12) 
fc=o "^^sT 



Since all the sets A™, /c = 0, . . . , m2™ — 1 are bounded (as subsets of -6(0, m + 1)), there exists 
for every m > 1 and every A;G {0, . . . , m2" - 1} a finite "firewall" (5f C M'^ (see [5] or Lemma 4.3 
in 131 and note that C (^rL™/2 := {xG M'^' : d{x, {Aff) > em/2}) such that ' 

Vn>l, VxgA^, (i(x,a„U/3r) = d(a;,(a„U/3r)nI^). 
Set = Uo<fc<m2--i/3r- Then, for every fee {0, . . . ,^2*" - 1}, for every xS 

d{x,an) > d{x,anUPr) = d{x,{anU ^T) > d(x, (a„ U Z?*") n I^). 

Set temporarily := card((a„ U Pi A"/^). First note that it is clear that 

nf card(a„ n Af) 



as n — > oo. 



n 



n 



It follows from the asymptotic L^-optimality of the sequence (q^) and the empirical measure 
theorem (see ()2.5() ) that 



card(a„ n AV') ^^'^ 
hm sup ^ — < — - — 2 

Jf~dXd 



n 



(2.13) 



so that 



n Jf'i+'-dXd ^ rn Jf^+rdXd 
hmmf — > ^ > — 2 



n n'i 



h^f—^dXd "^ + 1 V/^dAd 



On the other hand, for every k£ Im, 

(d(x, anWdXdix) > f d{x, {an U D n Af)dXd{x) > Xd{Af)eUAU{Af)) 



where U{A^) denotes the uniform distribution over ^™ (note that the inequality is trivial when 
Xd{A^) = 0). Then one may apply Zador's theorem which yields, combined with ()2.13p . 

s 

hminfnt / (d(x, a„))^dArf(x) > Ad(^n x liminf f x lim((n^)^e„„ ([/(^^)))^ 

( 



> UAf)x 



V 



m J f 'i+r dXd 

d 

/^m f^+^dXd 



m + 1 



J,,, X {Xd{At)y 



> J. 



> J. 



ni 



m + 1 



m 



m + 1 



fd+rdXd 



Xd{Al 



\ d fk + l 

f^+^dXd 



d+r 
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with the convention ^ = 0. 

Consequently, using ()2.12|) and the super- additivity of Hminf yield that, for every m > 1, 

s s 
( m \ d / f d \ d f ^ ^ s 

= Js,di^T) / f^dXA / fm{fm + 2-"'r — dXd. 



^m + lj \jRd J J{f>o} 

Now, by Fatou's Lemma, one concludes by letting m go to infinity that 

liminfni / id{x,an)YdP{x) > Js,d ( f f^dXa)' f f^~^dXd 



3 The upper estimate 



Let r,s £ (0,cxd), s > r. In this section we investigate the upper bound / d{x,anYdP{x) = 
0{n~^/'^) for L^'-optimal n-quantizers Un which is less precise than the lower bound given in the 
previous section. The reason for the restriction to (exactly) optimal n-quantizers instead of only 
asymptotically optimal n-quantizers will become clear soon. See e.g. the subsequent example 2. 
The formal restriction s > r is only motivated by the fact that the above upper bound is trivial 
when s <r since the LP{P) -norms are nondecreasing as a function of p so that 

s 

n'/'^ j d{x,anYdP{x) < [n''^'^ j d{x , OnY dP {x)y' . 

For a sequence (a„)„>i of finite codebooks in M!^ and b G (0, oo) we introduce the maximal 
function ip^ : M'^ ^ M+ U {oo} by 

Xd{B{x,bd{x,an))) . . 

'P»^^)-'^Pp(^B{x,bd{x,a^Y) ^ ^ 

(with the interpretation ^ := 0). Note that V'j, is Borel-measurable and depends on the underlying 
norm on M.'^. The theorem below provides a criterion based on these maximal functions that ensures 
the L^'-rate optimality of L^'-optimal n-quantizers. In Corollary ^ we derive a more applicable 
criterion which only involves the distribution P. 

Theorem 2 Let r,s £ (0, oo), s > r. Assume ^ and f \\x\\^~^'^dP{x) < +oo for some i] > 0. 
For every n > 1, let an be an L"^ -optimal n- quantizer for P. Assume that the maximal function 
associated with the sequence {an) satisfies 

for some b £ (0, 1/2). Then 

sup n'/'^ [ d{x, anYdP{x) < +oo. (3.3) 

n J 

Remarks. • In fact the proof of Theorem 21 provides a bit more information than ()3.3|) : it follows 
from (|3.6|) below and the (reverse) Fatou's lemma for sequences of dominated non-negative functions 
that there exists a real constant C3(6) > such that 

lim sup n 

'/'^ d{x,anY 

dP{x) < / lim sup n '/'^d{x,anYdP{x) < Csib) U^{xY^^'^^'^^dP{x) < +oo. 

n— >oo J J n— >oo J 
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• Condition ()3.2|) relates the upper estimate for / d{x, anYdP{x) to lower bounds for P-probabilities 
of (small) balls. If = P \i comes rather close to the minimal condition j f^^/('^+'^^dP < +00 
(see Theorem 1). Note that G L^{P) implies i/j'J ('^+''-) e L^{P) for every 6 > 6o- 

Proof. Let y G R'' and set 5 = 5n = d{y,an). For every x G B{y,6/2) and a G an, we have 

\\x — a\\ > \\y — a\\ — \\x — y\\ > 6/2 and hence 

d{x,an)>S/2> \\x-y\\, x £ B{y,6/2). 
Let /3 = /3„ = U {y}. Then 

d{x,P) = \\x-y\\, xG B{y,6/2). 
Consequently, for every b G (0, 1/2), 

CnAPy - en+iAPy > J d{x, anYdP{x) - J d{x, PYdP{x) 

> f {d{x,anY -d{x,PY)dPix) 

JB{y,5b) 

{d{x, OnY ~ \\x — y\\'^)dP{x) 

B{y,5b) 



> / {{d/2Y-{b6Y)dP{x) 

JB{y,5b) 

= {{1/2Y -bn6'P{B{y,b6)). 

One derives that 

^ P(W ^(^".'■(^)'' - ^n-,iAPY) (3.4) 

P[B{y,bd{y,an)) 

for every y G M'^, 6 G (0,l/2),n > 1, where C7(6) = {{1/2Y -U')-^ . Note that e„,r.(P)'"-en+i,r(i^)'" > 
for every n G N (see [3]). 

Now we estimate the increments en,r{PY ~ ^n+i,r{PY ■ (This extends a corresponding estimate 
in to distributions with possibly unbounded support.) 

Set en,r = en,r{P) for Convenience. Let {Va ■ a G a„+i} with Va = Va{an+i) be a Voronoi 
partition of MJ^ with respect to a^+i. Then P{Va) > for all a G a„+i and card a^+i = n + 1 (see 
0), 

a G an+i ■I \\x - a\\ dP[x) > — - } < 

I JVa n + 1 } 



card 

and 



card (a G a„+i : P(V;) > ^ 1 
I n + 1 J 

This implies that 

4e: 



< 



4 

n + 1 



:=|aGa„+i: / \\x - a\\'dP{x) < P{Va) < ^] 

I Jy„ n + 1 n + 1 J 



satisfies card/3„+i > (n + l)/2. Choose a closed hyper-cube = [— m,m]'^ such that Pr{K) > 3/4. 
The empirical measure theorem (see (|2.5|) above or [21 1^ for details) implies 

hm ^"^^(^^"^^ = PAR) 
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since Pr{dK) = Xd{dK) = 0. We deduce card(a„+i n i^) > 3(n + l)/4 and hence card(/3„+i n > 
(n + l)/4 for large enough n. Since one can find a tessellation of K into [(n + l)/8] V 1 cubes of 
diameter less than Cin~^/^, there exist ai,a2 G Pn+i, ^2 such that 

||ai -aall < Cm"^/'^ 

for every n > 3. Let 7 = a„+i \ {ai}. Using 

7) < ||x — a2|| < ||x — oill + II ai — a2||, 

one obtains 

enAPy - en+iAPy < J d{x, lYdPix) - J d{x, an+iYdP{x) 

= I \\x-aYdP{x)+ j d{x,-fYdP{x) - J2 I \\x-a\['dP{x) 

= f {d{x,jy -\\x-ai\\'')dPix) 

JVa^ 

< {r-l)f \\x-ai\rdP{x)+r\\ai-a2\\''P{Va,) 

^ 4(2'- - l)e;+,,, ^ 

~ n + 1 n + 1 

Consequently, using (|2.3j) . 

enAPy - en+iAPy < C^n-^'^^'^l^ (3.5) 
for every n G N and some finite constant C2 independent of n. Combining (|3.4|) and (|3.5() . we get 

VP(S(x,6(i(x,a„))/ 
< C3(6)^,(x)^/('^+^) (3.6) 

for every x G M'^,?! G N, 6 G (0,1/2) and some finite constant 6*3(6). The proof is completed by 
integrating both sides with respect to P . □ 

Next we observe that in case s < d + r a local version of condition ()3.2() is always satisfied. 

Lemma 1 Assume J Wx^ dP {x) < +00 for some r G (0,cxd). Let (a„) he a sequence of finite 
codebooks in M*^ satisfying J d{x,any dP{x) — > 0. Then the associated maximal functions ip,^ are 
locally in U'{P) for every p G (0, 1) i.e. 



VM,6g(0,oo), / ipPdP< +00 
Jb(o,m) 



B(0,M) 

Proof. Let M,b G (0, 00) and set A = supp(P). Then max^,g^(o,A/)nyl (^(^j c^n) ~* (see jS]) and 
hence 

C(M) := sup max d(x,a„) < +00. 

n>l x&B{0,M)nA 

One derives that 

B{x, bd{x, an)) C B{0, bC{M) + M) 
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for every x G 5(0, M) n A,n E M. Define the Hardy-Littlewood maximal function cp : —i- 
M+ U {00} with respect to the finite measures Xd{- D B{0, hC{M) + M)) and P by 

\d{B{x,p)f\B{Q,bC{M)+M)) 
P>o P{B{x,p)) 

Then 

Ai^) < fix) 

for every x £ B{0, M) n A. From the Besicovitch covering theorem follows the maximal inequality 

^ CiXdBiO,bC{M) + M)) 
P 

for every p > where the finite constant Ci only depends on d and the underlying norm. (See [H], 
Theorem|2ll9. The result in is stated for Euclidean norms but it obviously extends to arbitrary 
norms since any two norm on M*^ are equivalent.) Consequently, 



L 



B{oM)nA 



/roo roo 
ipPdP = / P{^P > t)dt < 1 + / P((/?P > t)dt 
Jo Ji 

/oo 
r^/Pdt < +00 

where C2 = CiXd{B{0, hC{M) + M). □ 

In case s > d + r under a mild assumption on the support of P a local version of (3.2) holds 
provided the minimal condition is satisfied locally. 

Lemma 2 Set A = supp{P). Assume P = = f-Xd, J \\x\\'''dP{x) < +00 for some r£ (0, 00), 

Xd{- n A) is absolutely continuous with respect to P and A is a finite union of closed convex sets. 
Let (an) be a sequence of finite codebooks in M'^ satisfying 

d{x,anTdP{x) 0. 

Then for p G (1,00], the associated maximal functions ipb are locally in LP{P) provided 1/f is 
locally in L^{P). 

Remark. The absolute continuity assumption Xd{- Ci A) <^ P does not follow from the absolute 
continuity of P: set P = clu-Xd where U = U„>o(r„ - 2-("+i), r„ + 2-("+i)), Q = {r„, n > 0} and 
c = l/Xd{U). Then supp(P) = M but ^ P since XdiU") = +00 and P([/^) = 0. 

Proof. Let M, 6 G (0, cxd) such that P(0, M) n ^ / 0. We have 

C = C(M) := sup max d(x,an) < +00 

n>l xeB{0,M)nA 

(see jSI) and hence 

(x) < sup ^^(^^^'^^^ 

for every x G P(0, M) H A. There exists a constant Ci G (0, 00) such that 

Xd{B{x, p) n P(0, M) r\A)> Cip'^ 
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for every x G 5(0, M) n A and p < bC since 5(0, M) n A is finite union of convex sets (see [3]). 
Define the maximal function if = tpf-i '■ M+ U {oo} by 

_ /b(.,p) f-'dP{.nB{o,M)) 
^^^^■"p>? p(i?(x,p)n5(o,M)) 

and note that dX^i- H A)/dP = 1//. One derives that 

A,(5(x,p)n5(o,M)n^) 
^^^"^ - P(5(x,,)n5(o,M)) 

for every x G 5(0,M)n A with C2 = Ad(5(0, 1))/C. By the LP(P(- n 5(0, M))-boundedness of the 
maximal operator ^ Vf-'^ (^^^ 13) Theorem |21 19), we obtain 

IIV'blB(0,M)||LP(P) < C'2||'/'1_B(0,M)||lp(P) < C'3||/~"'"1b(0,M)||lp(P)- 

This yields the assertion. □ 

Remark. One can replace the assumption on A = supp(P) by a local "peakless" assumption, 
namely 

-y , . ^ n • f Arf(5(x,^)n5(0,M)nsupp(P)) 
VM, c>0, mt TTTTT >0 

xeP(o,Af)nsupp(P),o<p<c X(i{B{x, p)) 

which can be satisfied by many subsets which are not finite unions of closed convex sets {e.g. if 
A = C^, C convex set with a non empty interior). 

An immediate consequence of Theorem |21 Lemma ^ and Lemma |^ concerns distributions P 
with compact support. 

Corollary 1 (Compactly supported distributions) Assume that supp{P) is compact and P'^ 7^ 0. 

(a) For s£ (r, d + r), the assertion \3. of Theorem\^ holds true. 

(b) Let s > d+r. Assume that supp{P) is a finite union of compact convex sets and P = P°' = f.Xd, 



J f-s/id+r)^p ^ 



Then the assertion of Theorem\^ holds true. If furthermore / > e > Xd-a.s. on supp(P), then 
the assertion i'j. 'A) of Theorem\^ holds for every s G (0, 00). 

Note that if s = d + r, claim (6) and the monotony of s 1— > || . ||s hold true if there is a positive 
real 5 > such that J f^^dP < +00. 

For distributions with unbounded support we provide the following condition. 

Corollary 2 (D^str^but^ons .Uh unbounded supports) Let s^ (r^d+r). Assume P'^ ^ 0, J ^r^dPi 
+00 for some 6 > and 

for some M,b £ (0,oo). Then the assertion il^.Aj) of Theorem\^ holds true. 
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Proof. Let xq £ supp(P). Then (i(xo,a„) — > (see 0). For ||x|| > N := \\xo\\ + sup„>i d(xo, a„), 
we have d{x, a„) < 2||x|| for every n > 1 and thus 

\d{B{x,t)) 

ip,{x) < sup 

i<26|lx'|| P{B{x,t)) 

for every x G B{0, NY,b > 0. The assertion follows from Theorem |^ and Lemma^ ^ 

For distributions with radial tails we obtain a condition which is very close to the minimal 
condition J /""Z^'^+'-UP < +oo. 

cou.. 3 H P . / < 

5 > and f = h{\\ ■ \\q) on B^^.^^^(0, with h : (i?, oo) nonincreasing for some R G and 

II • llo any norm on M'^. Let s G (r, d + r). // 

f{cx)-'/^'^+'UP{x) < +00 (3.8) 

/or some c > 1 ('and subsequently for every c' G (l,c]j, i/ien Assertion i'S. '■S\) of Theorem\E holds 
true. 

(6) Assume d = 1 and s G (1, 1 + r). Assume P = /.A, / |x|^"'"^dP(j;) < +cxd for some 6 > 0. If 
supp{P) C [Rq^co) for some Rq and f\{^R'^^oo) ^-^ nonincreasing for some R'q > Rq. Assume 

furthermore il5'.<^) /or some c > 1. Then the assertion o/ Theorem\^ holds true. 

Proof, (a) We may assume without loss of generality that || • || = || • ||o- For b G (0,1/2), let 
M = M{b) = N/{1 - 2b). Then for x G B{0,My,t< 2b\\x\\ and z G 5(0, t), we have 

\\x + z\\ > \\x\\ - \\z\\ > \\x\\{l - 2b) > N. 

Consequently, 



P{B{x,t)) = / f{z + x)dXd{z)= h{\\z + x\\)dXd{z) 

JB(0,t) JB{0,t) 

> h{\\x\\+t)Xd{B{0,t)) 

> h{\\x\\{l + 2b))Xd{B{0,t)) 

for every x G B{0,M)'', t < 2b\\x\\. Hence 

Xd{B{x,t)) ^ 1 ^ 1 

tim-W P{B{x,t)) - h{\\x\\{l + 2b)) /((l + 26)x) 

for every x G B{0,M)'^ and Corollary [2 gives the assertion. Item (b) follows similarly. □ 

Although Corollary |31 provides an efficient criterion for many families of multi-dimensional 
distributions since, in practise, most of them do have radial tails, it is also natural to provide a 
criterion which does not rely on this assumption. This is the aim of the next Corollary (which can 
also treat successfully most usual multi-dimensional distributions). 

c... . .™ P - / < _ . > „. „ 

that supp(P) has no peak i.e. 

inf UsnMP)nBix,p))^^ 

xGsupp(P),p>0 Xd{B{x,p)) 
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and that f satisfies the local growth control assumption: there exist real numbers e > 0, r/S (0, 
M, C > such that 

Vx, ye supp(P), ||x|| > M, ||y-x|| < 2r/||x|| =^ f (y) > C {f ix))'^+'^ . (3.10) 

f dP(x) ,__ 
Then, for every s G {r,jj^) such that / . < +oo (if any), the assertion i3. ^) from 

Theorem\^ holds true. In particular, if i!J.l(J\) holds either for e = or for every eG (0,e] {e>0), 
and if 

V.se (r,d + r), / dP{x)^ ^ j ^ ^^^^^^ 

JM.<i {f{x))'i+r J{f>0} 

then 123) holds for every s£ {r,d + r). 

Note that (if Ad(supp(P)) = +oo) Assumption ()3.9() is e.g. satisfied by any finite intersection 
of half-spaces, the typical example being M^^. Furthermore, a careful reading of the proof below 
shows that it can be slightly relaxed into: there exists a real c > such that 



r Ad(supp(P) n B{x , p)) 



'.= mf \ TTTTT vi ' ^ < P < c\\x\\ \ >0. 



Proof. Let xqG supp(P). Then, for large enough n, Un n B{xq, 1) / 0. Hence 

d{x, Un) < d{x, an n B{xq, 1)) < ||x — xoll + 1 < ||x|| + ||xo|| + 1 

so that d{x,an) < 2 ||x|| provided ||x|| > ||xo|| + 1. We will assume from now on (without loss of 
generality) that M > ||xo|| + 1 in Assumption ()3.1U|) . 

Then for every xG supp(P), ||x|| > M and every [0,2?7||x||] 

P{B{x,p)) = / /(2; + y)i{supp{p)}(a; + y)rfy 
Jb(o,p) 

> C(/(x))i+^Arf(i?(x,p)nsupp(P)) 

> CKf{f{x))'+'Xd{B{x,p)) 

where we used (|3.9|) . (|3.1()|) and ||x|| < 2r/||x||. Then 

Xd{B{x,p)) ^ 1 1 



P<2^\\ P{B{x,p)) ^ CKf{f{x)Y+^- 

Finally one concludes by Corollary [21 once noticed that 

fXd(B(x,p))\^ , 1 f dPix) 
sup ; ; '7x dPix) < — / -rnrr < +oo. □ 

^ll>A/p<2,,||x!| \P{B{x,p))J {CKf)—r J\\x\\>M (/(x))^ 

Remark. Note that the moment assumption (|3.11j) follows from the more natural moment as- 
sumption 

Vse(r, d + r), 3M>0,35>0 such that / \\x\\'^^^~^^ f{x)dx < +oo 

J\\x\\>M 

which is of course satisfied by all distributions having polynomial moment at any order. This follows 
from a standard application of Holder inequality (like in the remarks that follow Theorem ^ . 
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At the moment Gaussian distributions are most important for applications of the (r, s)-feature. 

Example 1 • Let P = N{0, S) be the d- dimensional normal distribution with positive defi- 
nite covariance matrix S. Then P satisfies condition (|ci.8|) from Corohary ^a) with h{y) = 
((27r)'^detS)~^/^e~^ ||a;||o := where || • || stands for the Euchdean norm on and 

M = 0. For s£ {r,d + r), choose cG (1, \/{d + r)/s). 

• For sS {d,d + r), the assumptions of Corollary |3Ja) are satisfied by the hyper- exponential distri- 
butions defined by 

/(x) oc exp(-a||x||^), xGW^,a,b>0. 



where || . || denotes any norm on (possible different of the underlying norm). Set h{u) tx 
exp(— au^). Note that if d = 1, then the normalizing positive real constant is given by = r{i/b) • 
In fact this holds true for all distributions of the form 

f{x) oc ||x||''exp(-a||a;||^), x G R'^, a, b > 0, c> -d, 

which e.g. include the (scalar) double Gamma distributions on the real line where 



f{x) = -r^bl^-^e-'^l^l, xgR, a,c>0. 

J \ J 2r(c) ' ' 

• As concerns scalar distributions let us first mention the Gamma distributions 

fix) = ^x''-^e-'^"l(o,oc)(x), X G M, a, 6 > 0, 

for which the assumptions of Corollary |31^6) are satisfied for every s€ (1, 1 + r). This holds as well 
for the Weibull distributions 

f{x) = bx'-^ exp(-x)^l(o,oo)(^), X € R, 6 > 0, 
the lognormal distributions 

f(^) = -7^exp(--^(logx - af)lm ^){x), xG M, a > 0, a G M. 

• The logistic distribution 



satisfies Corollary |21^a) and so do the symmetric p-stable distributions with p £ (0, 2) provided r < p 
ands< (j^(l+r). 

• For the Pareto distributions where 

/(x) = 6x-(^+i)l(i,^)(x), xgM, 6>0 

the assumptions Corollary 121^6) are satisfied provided that r < b and s < 6(1 + r)/{b + 1). 

• In fact all the above examples of distributions also fulfil the criterion proposed in Corollary ^ 
(although it turns out to be sometimes slightly more demanding). On the other hand some distri- 
butions need naturally to call upon this criterion, especially when their density is not monotonous. 
Let p, = g.Xi be an absolutely continuous probability distribution on [0, 1] satisfying 

< C_ < 5 < C+ < +0O. 
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Let {Yn)n>i denote an i.i.d. sequence of /i-distributed random variables and let N denote a Poisson 
random variable with parameter A > 1, independent of {Yn)n>i- Then set X = N + Y^r and 
P = P . One easily checks that 



= f\ with fix) = e-^-—^ g{x - (x). 



Let e > and rje (0, (1/2) A e) and xe R+, ye [x{l - r/),x(l + r/)]. As soon as n > [A] + 1, the 
sequence ^ is decreasing (to 0). Hence, if x > (1 + r/)[A] + 1 

fiy) > _C_^,,Al^l([x]!)i+^ 



(/(x))i+^ - (C+)i+^ [y]l A{i+^)N 



(C+)i+- (l + [x(l + r?)])!- 

One concludes using the Stirling formula that 

liminf inf , J^^}^^^ = 

On the other hand, Assumption (|3.11|) is fulfilled since X has finite moment at any polynomial order 
(see the remark following Corollary Note that when the density function g is not non-decreasing 
/ cannot be a non-decreasing function so that Corollary |31 does not apply. 

The exact optimality assumption made on the sequence (a„) in Theorem [21 (and Corollary^ 
is critical to get the optimal rate n~^^^. This is illustrated by the following counter-example. 

(Counter-) example 2. Let P = U{[0, 1]) and for n > 2 and -d G (0, oo) set 



= -=^2^}^{;^+(^-;^) '%n-l)' --' = '^ 



Let r G (0, cxd) and assume •& > -n^-p:- Using a non-Voronoi partition gives the upper estimate 



(l+r) 

E\X-X^-r < / l[o,„-.](x)|x-^rdx 



n „ 



^n^^^ n^' 2(n-l) 



dx 



r + l^2n^^ ^ nfi' (r -h l)2''(n - 1)^ 



Hence 



lim sup n''E| X - X"" | < lim sup 

n— >oo n— >oo 

1 



1 



2'"(r + i; 

so that in fact 



2'^(r + l)n('-+i)^-^' (r-hl)2^ 

■ = Jr^i = Qr{P) 



lim n'^E|X -X""r = Qr{P) 
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It follows that the sequence {an{'i!^))n is an asymptotically L''-optimal n-quantizer for every G 
(r/(r + 1), oo). Now, let s > r and G (r/(r + 1), s/(s + 1)). Then 

rl/2n^ 1 1 

so that 

n'ElX - X'^-l' > ——^ ^n^-'^Cs+i) 

' ' - 2«+i(s + 1) 

Consequently, 

lim n'ElX = oo. 

n^oo 

We finally comment on what can be said without condition (|H.2|) . 

Remark In the situation of Theorem [21 but without assuming ^^.'1^ let P = P'^ and s < d + r. One 
can deduce from ()3.6() and differentiation of measures that 

limsupni/'^d(- ,a„) < C3(6)^/"/~^/('^+'^) < +oo P-a.s. 

n— ►oo 

This improves considerably for absolutely continuous distributions and exactly L'^-optimal quan- 
tizers and a.s. result in Furthermore, by Lemma 1 and Fatou's lemma (also without (|3.2|) ). 

sup lim sup n'*'''^ / d{x,anYdP{x) < / limsupn*^'^(i(x, a„)''dP(x) 

M>0 n^oo Jb{0,M) J n-^oo 

< C^ib) f f~'/('^+'Up 



so that the minimal condition / f~^/('^+^')dP < +oo already gives the second inequality of the 
remark following Theorem [21 In particular, by Holder's inequality, the moment condition for P 
in Theorem [21 implies the minimal condition for s < r. The first inequality in this remark always 
holds for s = r since by Theorem 4.5 in 3 and (2.3), 



sup lim n''/'^ / d{x,anYdP{x) = sup Q,.(P)P^(5(0, M)) = Q,.(P) 
M>0"^°° Jb(o,m) M>0 

= lim n''/'^ / d(x,anYdP(x). 

n— >cxD J 



4 The critical and super-critical cases 

In this section, we investigate the upper bound for the rate of convergence when s > d + r. 

Proposition 1 Let r£ (0, oo). Assume 7^ and J \\x\\''''^^dP{x) < +00 for some 6 > 0. For 
every n>l, let an be an L"^ -optimal n-quantizer for P. 

(a) The critical case s = d + r: Assume there is a real > such that for every •&£ (0,?? A s), 
there exists a real number M > such that 

( sup WxfdPix) < +00. (4.12) 

B{0;Mr Vr<2fe||x|| P{B{x,r)) J 

Then, 

Vee (0,1 + r/(i), limsupn^+i"^ / (d(x, a„))^+''dP(x) < +00. 

n jRd. 
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(6) The super-critical case s > d + r: Assume there is a real ■& G {s - {d + r), s) and a real 
number M > such that 



Ib< 



I \ ^^^^^Y" \\xrdP{x)< +00. (4.13) 

lB{0;Mr \r<2b\\x\\ P{B{x,r)) J 

Then, 

supn^ / {d{x,an)ydP{x) < +00. (4.14) 

n JRd 

Proof. (6) It follows from 1)3. 6|) that, for every n > 1, 6g (0, 1/2), 

VxG (d(x,an))'^+" < Cbi',{x)n-^^+i\ 

On the other hand, xq £ supp(P) being fixed, for large enough n, a„ n B{xq, 1) 7^ so that 

yxe W^, d{x,an) < (1 + |x - xo|)l{|a;-xo|>i} + 21{|a;-x„|<i} < Co{\x\ V 1). 
Let tfj^ be the maximal function associated to (a„)„>i and b£ (0, 1/2). Consequently 

VxgM^ {dix,an)r < idix,an)r'''iCoi\x\Vl)f 

< Cln^ {i;Jx))^ (\\x\\ V if . 

It follows from Lemma ^ that, for every real number M > and 6 > 0, 

{ipf^{x)) \\x\\'^dP{x) < +00. 

B{0;M) 

3-1} 

Combined with Assumption ()4.13|) . this clearly implies V'if^' ^ 

(a) The above computations are still valid when s = d + r. So one concludes by considering 
arbitrarily close to which is precisely made possible by Assumption ()4.12|) . <^ 

From this Proposition one easily derives some corollaries in the formerly mentioned settings. We 
give them in details for the super-critical case. The adaptation to the critical case is straightforward. 

Corollary 5 Assume the global assumption on P in PropositionU\ holds and that P = f-Xd- Let 
t?G {s - {d + r),s). 

(a) If f is radial outside a compact subset o/M"^ (as defined in Corollary^ and if 

f{cx) d+r ll^ll dP{x) < +00 for some real c > 1 

then Assumption j^. is fulfilled. 

(b) If f satisfies Assumption ici. !A) and Assumption \3.1U\) for some parameter e > and if 



/ f{x) '^+'- ||x|| (iP(x) < +00 



then Assumption i4.13{) is fulfilled. 



The above assumptions are clearly fulfilled by the normal distributions (and in fact most dis- 
tributions mentioned in Example 1). Numerical experiments seem to suggest that the critical rate 
n^~^d cannot be improved for P with unbounded support, that is. 



Vs>0, limn^+d / d(x,anYdP(x) = +00. 
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5 A sharp rate for distributions on compact intervals 



Proposition 2 Let P = f.\^ where f : [a, b] — > M+ is a Lipschitz continuous probability density 
function, bounded away from on [a,b]. Let r > 0. Let (an)n>i be a sequence of L^ -stationary and 
asymptotically L^ -optimal [a, b]-valued n-quantizers. Then for every s > 0, 



lim ^n"* j mill \x - aYdP{x)^ = Qr,s{P)- 



Note that we do not assume a priori that the n-quantizers are L''-optimal but only stationary 
and asymptotically L^-optimal. 

Proof. It is possible in this scalar setting to number the elements of a n-quantizer with respect to 
the natural order on the real line. Furthermore, one may assume for large enough n that an C (a, b) 
and carda„ = n. Namely, we set 

On = {a^, . . . , a^}, a < < < • • • < a^"-^ < a"^ < b. 

We also set = a and a^"*"^ = b for convenience. For every n > 1, set for every xG [a,b] 

n 

ipr,n{x) = n ^ Aa^ly^(„„) (x) (5.1) 

k=l 

where Vk{oin) denotes the Voronoi cell of and Aa^ = — a^"^, A;G {l,...,n + l}. We know 
from j|4| (see the proof of Theorem 4(a), p. 1101) that h being Lipschitz continuous and bounded 
away from on [a, 6], there is a positive real constant Cafij such that 

max min Aq„. (5.2) 

l<fc<n+l ~ ' '■' l<fc<n+l " ^ ' 

1 " 

A proof of this last inequality can be found in the appendix. Set P^.n = — ^ the empirical 

^ k=l 

measure associated to the n-quantizer an- By the empirical measure theorem (recalled in 1)2. 5() ). it 
follows that for large enough n, 

\/a,b'e[a,b], — i— L_LL — > Pri[a ,b']) := Cr ^ f^'+^dX^ 

n -I'T+i J^i 

uniformly with respect to {a' ,b')^ [oj^]^ (where Cg^s = {f^ g"^ dXi)^^). Combining this convergence 
with (|5.2|) and using that / is bounded away from on [a, b] implies that there is a real constant 
J, J such that, for every n > 1 and every a', 6' G [a, b] 

h' — ' C 

< min Aa^ < max Aa^ < ^^(6' - a'). 

C'abf''^ l<fc<n+l,ofee[a',fe'] l<A:<n+l,a^G[a',fe'] n 

A straightforward application of the Arzela-Ascoli theorem yields that the sequence {(pr,n)n>i is 
relatively compact for the topology of uniform convergence and that all its limiting functions are 
continuous. Let ipr denote such a limiting function. On the one hand, for every a', b' £ [a,b] , 



nb' " 

J a' 1 



— a 



,/] j — > u — a as n ^ 00. 

fc=i 
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i.e. (pr,n-Pr,n converges weakly to the Lebesgue measure over [a, b]. On the other hand (pr,n-Pr,n 

converges weakly to (pr.Pr = Cr_i_(prf''+^.X-^ since tpr is continuous. Consequently, ipr is entirely 
and uniquely determined by this equation which implies that 

u 



Vr,n ^ y^r ■= \ / f^'+^dXi ) / ''+1 aS n ^ GO. 

Now for every A; = 2, . . . , n — 1, 

/ \al-xYf{x)dx-f{ai)j \al-xYdx < n^[/U / W^-xY+^dx 

JVk{a^) JVkia^) JVk(ah) 

= [/] 



n 



n 
s + 2 



+ 



[/] 



Lip 



(s + 2)2^+2n2 



2 J \ 2 

s+2 



and 



n 



x\^dx = 



is 

Jl,s 

'2n 



since 



[/] 



Lip , 



n 



For k = 1 and n, the above equalities are no longer true due to edge effects but both induced 
errors are 0(l/n). Consequently 



n 



E\X -X^ 



2 n 



fe=l ^ 



n 



On the other hand 



2 / {ifrY'^^fdPr as ra ^ GO. 



Finally 



[\iPrr+^fdPr = ( f'f^^dX,] f' f'-^^f^^dX, = ( f'f^^dX,] f'f-^^dX, 
J a \Ja J J a \Ja J J a 
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which completes the proof. ^ 

Remark. In fact, we proved a shghtly more precise statement, namely that, for every r,s > 0, 
E\X-X''-\' = ^^:^ + o(^) as n^oo. 

6 Application to quadrature formulas for numerical integration. 

6.1 Numerical integration on 

Thanks to this extension, it is now possible to take advantage of the stationarity of quadratic 
quantizers for a wider class of functions. 

Let / : M'^ ^ R be a twice differentiable function such that 

\\D^f{x)\\<A{\\x\\''~^ + l) 

for some non-negative real constants A and B and r]£ (0, d]. Let X be a r.v. with distribution P = 
W^. Assume that Xs L'^~^'^{W) and that for this distribution the extended exact L*-quantization 
rate of a sequence of optimal quadratic n-quantizers (an)n>i holds true for every s £ (0,2 + d). 
First, let a be an optimal quadratic quantizer and X" a related Voronoi quantization of X. For 
every couple {p, q) S (1, +oo)^ of Holder conjugate exponents, one has 

|E/(X) -E/(X")| < ^||dV(^")IUI \\X-Xm- (6.1) 

Note that this error bound follows from the stationarity property of a quadratic quantization 
that is E{X \ X'^) = X'^ since 

|E/(X)-E/(X")| = |E/(X)-E/(X")-E(Z)/(X").(X-X"))| < ^E{\\D^f{X")\\ ). 
Now, set p = ^ and g = |±| in Equation Then {d-r])p = d + 2 and 2q < d + 2so that, 



-oo 



and 

||X = 0{n-i). 

Consequently the optimal rate 



|E/(X)-E(/(X"-))| =0(n 



for numerical integration by (optimal quadratic) quantization holds for a wide class of (twice differ- 
entiable) functions whose growth at infinity can be infinitely faster than quadratic (the quadratic 
case is obtained for q = 1 and p = +oo for functions having a bounded Hessian; then there is no 
need for the extended quantization rate). 
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6.2 Numerical integration on the Wiener space 

Let 



n>l 



the Karhunen-Loeve {K-L) expansion of a standard Brownian motion (VFt)tg[o,T]- (e^,A„)n>i is 
the Karhunen-Loeve (orthonormal) eigensystem of the (nonnegative trace) covariance operator of 
the Brownian motion C^(s,t) = s A t, indexed so that A„ decreases (to 0). All these quantities 
do have explicit expressions (see [Him among others). In that K-L expansion, (^n)n>i is then 
a sequence of i.i.d. random A^(0; l)-distributed random variables. One considers a sequence of 
(scalar) product quantizations of the Brownian motion W i.e. 



W 



fc=i 



where for every iV > 1, iVi, . . . , A^fc > 1, 1 < A; < and TVi x • • • x < N, is an L^-optimal 
iVfc-quantization of (see [71 |H] for more details). Then denoting | . 1^2 the quadratic norm on 

L'^{[0,T],dt) one has for every [2,3), 



\W-W 



e(\w-w 

(E{\W-W\l2 
\\\W -W\l2 



T 
1 
2 
s 

T ^ 



k=l 



< 



+ E 



,fc=i 



E^fcii^fc-4^ii'+ E ^'^ii^' 



,fe=i 




where we used that | > 1 so that 



J, is a norm. Consequently, using that the normal distribution 
satisfies Theorem |21 (see Example 1), for every s£ [2,3) there exists a real constant Cs such that 

ir-^^^^ii2<iiH^-w^^iL<cJE^+ E =CsW-y^''\\.- (6.2) 

\k=\ k k>m^+l J 

Consequently, if furthermore {W^)n>i is a rate optimal sequence for quadratic quantization, then 



\W-W 



\w-w 



Ciy(logA^)" 



(6.3) 



Note that the result clearly holds for s G (0,2] (simply using the monotony of the L*-norm). 
Furthermore, except for the final rate 0((log A^)~2) for which we refer to 0, the Brownian motion 
plays no specific role among Gaussian processes : the inequalities in (|6.2() extend to the product 
quantization of any Gaussian process iXt)t£[o,T]- 
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In fact, as emphasized in IS], optimal scalar product quantization is rate optimal but it is also 
possible to produce some rate optimal sequences based on d-dimensional block product quanti- 
zations. In fact it is shown in |Hj that as d grows one can obtain the sharp convergence rate of 
high-resolution quantization (it is even possible to choose this dimension d = d{N) as a function of 
to achieve this sharp rate) . From a numerical viewpoint producing these d-dimensional marginal 
blocks is more demanding than producing scalar optimal quantizations, but some issues are in 
progress on that topic. 

The point here is that, if one considers some d-dimensional marginal blocks to produce a rate 
optimal sequence of quantizations of the Brownian motion, such a sequence will satisfy ()6.3|) for 
every s€ (0, d -|- 2) instead of (0, 3). This seems to be an interesting feature of d-dimensional block 
quantization. 

The practical application to numerical integration on the Wiener space is formally similar to that 
on MJ^ except that one considers some functionals F on (L^([0, T], dt), | . 1^2 ) instead of functions / 

defined on M"^. Then, one approximates E F{W) by E F{W^) having in mind that a closed formula 
exists as well for the distribution of {the weights, see jTT]) which makes the computation of 
KF{W'^) fully tractable (with negligible cost). Then, if the Hessian D^F of F satisfies 

yuj€ L'^{[0,T],dt), \\D^F{uj)\\ < ^(|cj|22 +1) 

for some positive real constant A, r > similar computations as those carried out in the d- 
dimensional case implies that 

\EF(W) -EF(W^)\ < 

logN 

Without the extended quantization rate 1)6. 3|) such a rate only holds for functionals F with a 
bounded Hessian. For further applications (like Romberg extrapolation) and numerical experiments 
devoted to path-dependent option pricing, we refer to [TT] . 



Appendix 

In this appendix we closely follow the proof of Theorem 4 in ^1. This part of the proof is 
reproduced for the reader's convenience. 

Set Aa^ = — a^~^. After an appropriate change of variable the L^'-stationarity Equa- 
tion (see 13]) 

'^ Lk+^h-i - xr"^sign(a^ - x)/(x)(ix = 0, 1 < k < n, 

•' 2 

reads: 

J a4 ^''"^sign(?;) /(a^ + v) dv = 0, 1 < k < n. 
Introducing /(a,^) we obtain: 



2 



r2^ 



Aq 



fe+1 



v^-' {f{ai + v)- f{ai)) dv + /(a^)L_n_^. 
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Setting V := and v := — u respectively we have, since f{a^) > 0, 
(Aaj(l + ry^n —^^ du 



■f 
Jo 



{Aat'Y I 1 + r / — J^^^ du 



Setting y) := / ru^ ^— — ^ — 77^ ^ du finally leads to 

Let [a, b] C (m, M) and let L^''' denote the Lipschitz coefficient of / on [a, b]. 

v^, e'G [a,b], \Hic,a\ < yJ^'f < c^a,6,/ le-e'i 



(6.4) 



since / is bounded away from on [a, 6]. Consequently: 

\H{atai~^)\ < Ca,bj\ai " «n"^l and |/7(a^, a^+^)| < C„,6,/|a^ - Q^+^|, 
whenever lie in [a, 6]. 

The p.d.f. f being bounded away from on [a, b], inaX|^ ^ a<a''='=^<b} ^^^(Aa^, Aa^+^) goes to zero 
(see e.g. Proposition ?? in [1]), so we can estimate the right hand of (|6.4|) : 

= ^MHiatat') " + 0(max(Aa^, Aa^i)^) 

1 + -H (a^, On ) 

where 

max |0(max(Aa^,Aa^+i)2)| < C max max(Aa^, Aa^+i)^ 0. (6.5) 

{k I a<OLi^^ <h} {k I a<a^^ <h} 

Now, we have ^ Aa^ < b — a, so that: 

{fc / a<a'^'^'^<b} 

^ (Aa^)2 < (5-a) sup (Aa^) = (6 - a)o(l). (6.6) 

{fc/a<a;=±i<fe} {fc/a<aS±i<fe} 

Now set /i£ := exp(i/(a^+^ a^+^+i)-i/(a^+^ a^+^-i)) and := 0(max(Aa^+^+i, Aa^+^)0. Note 
that, as long as and af^^ G [fi,^], miiii<£<p-i /i^ is bounded away from by a real constant 
h>0 not depending on n since the function H is bounded away from —00 on [a, 6]^. Hence 

^ <n(^^+i<5^i) = (nMn(i+i^)<(nMexp(ET)^(n^^)e^'-'^^°^'^- 



=1 



One obtains a lower bound the same way round. Assume that n is large enough so that, maxjj^, a < 
a^rt^ <b} < h/2. Using that ln(l - u) > -2u if < u < 1/2, 

f^TT > ( n n (1 - f ) > ( n -P(-2 E f ) > ( n /^^) ^^^-'^^ • 
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Thus it yields that for large enough n and every a'^_^_p, a" E [a, b]: 

exp(^ Hhi) -{b- a)\o{l)\) < < eMT. Hhi) + {b - a)\o{l)\). 

1=1 V^«" / e=i 

On the other hand, using once again that / is bounded away from on [a, 6], yields 



p-i 

E Hht 

1=1 



< L 



a,b r 



r + 



p-i 
£=1 



< C(6-a). 



Combined with Inequality (|6.7() . this provides 



(6.7) 



(6i 



max Aa„ < C min Aa„. 

{k / a<x'i'^\<a*i<b} {k / a<a^-^<a^<b} 
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